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HQET lagrangian up to 1/m^ terms is discussed. Consequences of reparameterization invariance 
are considered. Results for the chromomagnetic interaction coefficient at two loops, and in all 
orders in the large-/3i approximation, are presented. 



1 HQET lagrangian 

. QCD problems with a single heavy quark staying approximately at rest can be conve- 
Q\ \ niently treated in the heavy quark effective field theory (HQET) (see for review and 
references). We shift the energy zero level: E = m + uj, and consider the region where 
^ ' residual energies uj and momenta p are not large: ~ |p| ~ A ^ m. The effective field 
^ ■ theory is constructed to reproduce QCD on-shell scattering amplitudes expanded to some 
order (A/m)". This is achieved by writing down the most general effective Lagrangian 
consistent with the required symmetries, and tuning the coefficients to reproduce QCD 
on-shell amplitudes. Terms with D^Q can be eliminated by field redefinitions. 
The most general lagrangian up to is 0-0 



> 

Cj ; L = Q+iD^Q 

+£Lq+d'Q + ^Q^B . aQ + ^Q^{D x E - E x D) ■ aQ + ^Q^[D ■ E]Q 
2m 2m 8m^ ^m^ 

0\ 



+ ^Q+D^Q + B-a}Q- ^^Q^D^B • BD^Q (1) 



Dh: +^g+(^fi ■D + D-BD)-aQ + P^g+(^ ■[DxB] + [DxB]- D)Q 

■ - E')Q - ^^Q^E'Q + Tr(52 - E')Q - -§^g+ Tr E'Q 

>< : +^g+(5 xB-ExE)-aQ- p^^Q+{E x E) ■ aQ + ■ ■ ■ 

where Q is 2-component heavy-quark field. Here heavy-light contact interactions are 
omitted, as well as operators involving only light fields. 

HQET can be rewritten in relativistic notations. Momenta of all states are decomposed 
as p = mv + k where residual momenta A; ~ A. The heavy-quark field is now Dirac spinor 
obeying ^^Q^ = Q^. The lagrangian is 



= Qjv ■ DQ, - -±Q^DlQ, - -^Q^G ,,,0^'' (2) 
2m Am 



+^Qv{Dl, G^^jv^a.^Q. - ^Qy^[DlG,,]Q. + ■■■ 
8m^ 8m^ 

where D± = D — v{vD). The velocity v may be changed by an amount Sv < A/m without 
spoiling the applicability of HQET and changing its predictions. This reparameterization 
invariance relates coefficients of varying degrees in 1/m • 

At the tree level, there are easier ways to find the coefficients Gi than QCD/HQET 
matching: Foldy-Wouthuysen transformation |14, 15|, or using equations of motion |^ 
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(or integrating out lower components [0, |T^) followed by a field redefinition. The result 
is 



Ck — Ci 

Cw2 



Cpip 



Cd — Cs — Ck2 — Cwl — Cal — Cbl — 1 , 
= Cm = Ca2 = Ca3 = Ca4 = Cf)2 = . 



(3) 



However, these algebraic methods don't generalize to higher loops. 

At 1/m level, the kinetic coefficient = 1 due to the reparameterization invariance 0. 
One-loop matching for the chromomagnetic coefficient Cm was done in 0; two-loop 
anomalous dimension of the chromomagnetic operator in HQET was obtained in [|18[ [T^ , 
and two-loop matching was done in ; in , all orders of perturbation theory for Cm 
were summed at large Pi. 

At level, the spin-orbit coefficient Cg = 2Cm — 1 due to the reparameterization 

invariance |2T|-|^. The Darwin term reduces to a contact interaction. One-loop match- 
ing for the heavy-light contact interactions was done in [24|. The one-loop anomalous 



dimension matrix of dimension 6 terms in the HQET lagrangian was obtained in |T5 



22 - 25 



At level, one-loop matching was done in ^ for the terms involving the heavy-quark 
fields twice and the gluon field once. The one-loop renormalization of dimension 7 terms 
in the HQET lagrangian was recently considered [EH]. 



2 Matching quark— quark vertex 

Renormalized QCD on-shell quark-quark proper vertex 

— — m)u (4) 

gets no correction in the on-shell renormalization scheme. QCD spinors are related to 
HQET spinors by the Foldy-Wouthuysen transformation 

/UP \ 

Expressing QCD proper vertex via HQET spinors, we obtain 

Uy^u^ H (6) 

2m 

Let's denote the sum of bare 1-particle-irreducible self-energy diagrams of the heavy 
quark in HQET at 1/mP as — i^-j^S(a;), u = kv. At the 1/m level, self-energy diagrams 
with a single chromomagnetic vertex vanish. Let the sum of bare diagrams with a single 
kinetic vertex be —i^^^T.k{uj, k'^). Consider variation of at v v + 5v for an 
infinitesimal 5v (v 5v = 0). All factors can be combined into a single one, and the 
variation 6-/) in it provides the variation of the 7-matrix structure in front of S. There 
are two sources of the variation of S. Terms from the expansion of denominators of 
the propagators produce insertions ikSv. Terms from the vertices produce igf^dv^. Now 
consider variation of at — > /c_l + 5k ±_ for an infinitesimal 5k j^. Quark-quark kinetic 
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vertices produce i^kSk±; quark-quark-gluon kinetic vertices produce i^gt°-Sk'^; two- 
gluon vertices produce nothing. Therefore, 



dk'f 



9S 
) 



(7) 



into account |pr 



and 



This is the Ward identity of the reparameterization invariance first derived in |TD]. Taking 

^/c^, we obtain 



dk^, 



duj 



(8) 



The right-hand side does not depend on k'j_, and hence 



Sfe(u;, k] 



k]_ + Sfco(u;) 



(9) 



This result can also be understood in a more direct way. Only diagrams with a quark- 
quark kinetic vertex contain /cj^; its coefficient is is i 2;^. The sum of diagrams with a unit 
insertion is —i^- Note that diagrams with a quark-quark-gluon kinetic vertex vanish 
because there is no preferred transverse direction. 

On the mass shell (cj = 0), the renormalized HQET quark-quark proper vertex is ^Zqu^ 
[— fc5^ + Sfc(0, fc^)]'Ut, = —^Zq [1 — k'juyUy. On the mass shell, only diagrams with 

finite-mass particles in loops contribute (e.g., c-quark loops in 6-quark HQET) (Fig. |l]). 
Taking into account Zq^ = 1 — and comparing with (|^), we finally obtain 



Cfc(/x) = l. 



(10) 



This argument works for an arbitrary hence, the anomalous dimension of the kinetic- 
energy operator in HQET vanishes exactly. In a similar way, it is not difficult to prove 
that 

Ck2=l- (11) 
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Figure 1: HQET quark-quark proper vertex on the mass shell 
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3 Matching quark— quark— gluon vertex 



QCD on-shell proper vertex is characterized by 2 form factors: 



2m 



4m 



u{p) , 

,2 



(12) 



e{q^) = l+e'^ + ---, fi{q^) = fi + fi' 



+ 



m^ 



m^ 



The total colour charge of a quark £:(0) = 1 due to the gauge invariance. Ward identities 
in the background field formalism ||2^ are shown in Fig. where the large dot means 
convolution with the gluon incoming momentum q and colour polarization e"", the second 
equalities are valid only for an infinitesimal q (or in the case of an abelian external 
field), and {t"')^'^ = if'^'^^ in the adjoint representation. Therefore, the QCD proper vertex 
^1l{Pi1) = obeys A^g^e'^ = + ge"t") + S(p) for infinitesimal g, or A^(p, 0) = 

— ^g^y . The form factor is projected out by e{Qi) = Zq[1 + ^ Tr A^t>^(l +|^)]. On the mass 
shell, I Tr = (1 - Zq^)v^, and hence £(0) = 1. 



n 
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+ 
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Figure 2: Ward identities in the background field formalism 



Let's denote the sum of bare vertex diagrams in HQET at 1/mO as i^tVi±^[l + A(u;, A)], 
where A = qv = lv' — uj. The Ward identity for the static quark propagator is the same 
as for the ordinary one (Fig. |^). Therefore, Ae°'t"'A{uj, A) = — S(u; + Ae°'t"') + S(u;) for 
infinitesimal A, or 



A(tu,0) 



duj 



(13) 

^(-+^)-^(-) exactly. The 



It is interesting, that for an abelian external field A(a;, A) — ^ 
total colour charge of a static quark Zq[1 + A(0, 0)] = 1, as expected. 
The 1/m HQET bare proper vertex has the form 

'S^^"^ [(1 + Afc)(P + P')! + + AkiPl + KiP'l + ^k2ql)v''] 
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+ »^9t"ii^l7",rflii^(l + A™), (14) 

where all Aj depend on A; A'^]^(u;, A) = Afci(u; + A, —A); Afc(a;, A) = Kk{uj + A, —A), 
and similarly for A^o, ^k2- Similarly to the previous Section, we can see that variation of 
the leading vertex function dX v v + 5v coincides with that of the kinetic-energy vertex 
function at p± — > p± + if 5v = ^6p±. This requires 

A,(a;,A) = A(^,A), A',,(^, A) = (15) 

(and hence Aki{uJ, A) = — ^) A(uj, A)). The Ward identities of Fig. |^ result in 

Ai,o{uj,0) = -^^^^ , A,,2(a;,0) = (16) 

(in an abelian external field, AkoiuJ,A) = _ s.o(^+A)-SfcoH ^ Afc2(cu, A) = 0). 
Reparameterization invariance relates the spin-orbit vertex function to the chromomag- 
netic one, but we shall not discuss details here. 
The on-shell HQET vertex at the tree level is 

Mk') ( + C,^^±^ + C^Mlll + c,f-^v^ + c}^v^ + ■■■') u^k) . (17) 
2m 4m 8m^ 8m^ 



As we have demonstrated above, there are no corrections to the first two terms. Other 
terms have corrections starting from two loops, if there is a finite-mass flavour (such as c 
in 6-quark HQET). Expressing the on-shell QCD vertex via HQET spinors, we obtain 



Uv{k') 



Therefore, the coefficients in the HQET lagrangian are 

Cfc = l, Cm = ^l, Cd = 8£' + 2/i-l, = 2/2-1. (19) 

The first one has no corrections (0). The coefficients ([T9|) are not independent: 

= 2C„ - 1 . (20) 

Probably, reparameterization-invariance Ward identities yield relations among corrections 
from finite-mass loops in HQET which ensure the absence of corrections to (^). However, 
we shall not trace details here. 

Similarly, at the 1/m^ level, the coefficients in the HQET lagrangian are 

C^i = V + i/i + |, C^2=V+i/i-|, Cp.p = /i-l, Ca/ = -4£'-i/i + i. (21) 

They are not independent: 

Cw2 = Ciui ~ 1 ) Cpip = Cm — 1 , Cm = I {Cm — Cd) . (22) 

Calculation of Ca, Cb requires matching amplitudes with two gluons. Calculation of 
contact terms requires matching amplitudes with light quarks. 
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4 Chromomagnetic interaction at two loops 



As we know, the kinetic coefficient Ckin) = 1, and the only coefficient in the HQET 
lagrangian up to 1/m level which is not known exactly is the chromomagnetic coefficient 
Kn(/i). It is natural to find it from QCD/HQET matching at /i ~ m where no large 
logarithms appear. Renormalization group can be used to obtain Cm at /i <^ m: 



CmifJ^) = Cm{m)exp 



where C^(m) = I + + C, {f^Y + 



( \ 

7m (a) da 



\ as{m) 



(23) 



In 



d log Zm 

dlog n 



7ii^ + 72 (^)' + ---is the 



anomalous dimension of the chromomagnetic operator in HQET, and the /3-function is 
= = /5i t + (t)' + ■ ■ ■ (where A = '^Ca - ITfUj). If L = logm/fi is not 

very large, it is better to retain all two-loop terms and neglect higher loops: 



CM = l + iC,-^,L) 



47r 



+ [C2 - (Ci7i + 12)L + 71 (71 - /9i) L^] (^)' • (24) 



This approximation holds up to relatively large L because C2 is numerically large. If L is 
parametrically large, then it is better to sum leading and subleading logarithms: 



Cm(/i) 



as{m) 



71 
2/3i 



1 + Ci 



An 



2/3? 



(25) 



In this case, we cannot utilize C2 without knowing 73. In general, the solution of (|23|) can 
be written as 



CM =CmK{fi), Cm = as{m)^^i{l + Sc) , Sc 



Ci^^ + C2 ' ^ ^ 



47r 



47r 



+ ■ 



(26) 



where Cm is scale- and scheme-independent. 
As a simple application, we consider B-B* mass splitting p8 

1 



1 



TUB* - ms 



2CmifJ') 2 



3m 



+ ^ [<^m(^)pL(/^) + Cm{p)pmmif^) " Cs{fi)pl{fi)] , (27) 



where Pmif^) Plif^) local matrix elements of chromomagnetic interaction and spin- 
orbit one, while plmip) and Pmmip) kinetic-chromomagnetic and chromomagnetic- 
chromomagnetic bilocal matrix elements (in the later case, there are two 7-matrix struc- 
tures, 1 and a^^; the coefficient of the second one is implied here). Introducing renormal- 
ization group invariants 



K{p^)fiUp) , Plm = K{p)plm{p) + [1 - K{p)\ pl{p) 



K\p)p: 



.3 A3 
mm ' rs 



{2i 



we can rewrite it as 



TUB* - itlb 



2C'm ^ 2 I 
'Pm "T 



3m 



3m^ 



^■m {Pkm ~ 2Ps) + C^Pmm + P'. 



(29) 
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In order to obtain Cm, we should calculate the heavy-quark chromomagnetic moment 
(Fig. P). All on-shell massive integrals can be reduced to 3 basis ones 



-'o 






(30) 



using integration by parts |3^-|3^- Iq ^^"^ -^i expressed via F-functions of d; I2 is 

(31) 



expressed via /g, Ji, and one difficult convergent integral p2 



/ = 7r'log2- -C(3) + 0(£). 



The result has the structure 



4 — 4e 

9om *%„2 



(32) 



H — (^^j^Y i^F^ CpCa, Cj^, CpTpTii, CaTfUi, CpTp, CaTf) x (Jq , Ji, I2) . 

Now we express it via as{fi) and expand in e. The coefficient of 1/e gives the anomalous 
dimension 

= 2Ca^ + ^Ca (UCa - ISTpUf) {^y + --- (33) 
The chromomagnetic interaction coefficient at = m is 



C„(m) = 1 + 2{Cf + Ca) 



Air 



+ 



20 



269 \ 



/ + - 31 + CpCa 77^ + ttVt' + — + 77^ - + ^ 



9 ) 



17 



805 



9 



27 



+CFTFni 



9 / V 9 27 



(34) 



T I ^6^2 ^^^^^ T 298 



47r 



1 + 1^^ + (21.79 -1.91n.)f 

6 TT V 



7r 



The coefficient of (as/vr)^ is about 11 for n/ = 4 light flavours. It is 40% less than the 



expectation based on naive nonabelianization |Q. The contribution of the heavy quark 
loop to this coefficient is merely —0.1. 



5 Chromomagnetic interaction at higher loops 



Perturbation series for Cm can be rewritten via f3i instead of Uf-. 

00 L— 1 



C„(/i) = 1 + V V aLnPia^ 



L=l n=0 



l + -^/(Aas) + 
Pi 



Mn H], is missing; in 



the leading logarithmic running of Cm{fi) has a wrong sign. 



(35) 
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There is no sensible limit of QCD in which jSi may be considered a large parameter 
(except, may be, nj — oo). However, retaining only the leading (3i terms often gives 
a good approximation to exact multi-loop results ||33|. This limit is believed to provide 
information about summability of perturbation series MM. At the first order in 
multiplicative renormalization amounts to subtraction of terms; 



(47r)2 ^ ! + /?/£' - 47r 21og/i/AMs 
The perturbation series (B^) can be rewritten as 



(subtractions) + O [ 

Pi 



Knowledge of the function F{e, u) allows one to obtain the anomalous dimension 



7r, 



21 



F(-/3,0) + O 



PI 



and the finite term 





CM = 1 + 



1 /,/(M)-F(0.0) 1 

Pi J e Pi 

-/3 



du e 



.,/^ F(0,.x)-F(0,0) 
u 



PI 



(36) 



(37) 



(3^ 



(39) 



(this method was used in |3^; see references in this paper). Renormalization group 
invariant is 



5c 



1 



due '3i"="5(n) + . 
Pi Jo VPi 



5(n) 



e 3 



, F(0,n)-F(0,0) 



(here is taken at = m in the V^-scheme, exp( 



U 



U 



■ Kv \ —2u\ 



(40) 



li=m 





LI I 




L 2 I 



b c 

Figure 4: L-loop diagrams with the maximum number of quark loops. 

The function F{e,u) is determined by the coefficient of the highest degree of Uf in the 
L-loop term, which is given by the diagrams in Fig. ^ Calculating them, we obtain 

/i\2« r(i + M)r(i-2M) 



F{e,u) 



m/ 



r(3 - u - e) 



-D{eY''-^N{e,u) 



D{e) = 6e^^r(l + e)B{2 -e,2-e) = 1 + U + 



(41) 



N{e, u) = Cf^u{1 + u- 2eu) + 



2-u-e 
9 



(2 + 3n - 5£ - Qeu + 2e^ + Ae^u) . 



This gives the anomalous dimension 



7r, 



C, 



as /j(l + 2/3)r(5 + 2/3) 
-271 24(l + /3)r3(2 + /3)r(l-^) 



(42) 



a. 



1 + 



13 A«£ 

6 Att 



+ 



This perturbation series is convergent with the radius 
r(M)r(l - 2u) 



Siu) 



r(3 



u 



[Au{l + u)Cf + i(2 - u)(2 + 3u)Ca\ - e 



< Att. The Borel image of 6c 



u 



(43) 



has infrared renormalon poles at u = ^. They produce ambiguities in the sum of the 
perturbation series for 5c, which are of order of the residues ~ (Ay/m)". The leading 
ambiguity {u = |) is 

8Cf) m ' 



(44) 



where Am is the ambiguity of the heavy-quark pole mass MB, EB 



Physical quantities, such as the mass splitting (|27|) , are factorized into short-distance co- 
efficients and long-distance hadronic matrix elements. In regularization schemes without 
a hard momentum cut-off, such as MS, Wilson coefficients also contain large-distance 
contributions which produce infrared renormalon ambiguities. Likewise, hadronic matrix 
elements contain small-distance contributions which produce ultraviolet renormalon am- 
biguities. In other words, the separation into short- and long-distance contributions is 
ambiguous; only when they are combined to form a physical quantity, an unambiguous 
result is obtained. Cancellations between infrared and ultraviolet renormalon ambiguities 
in HQET were traced in [R^ . 



Ultraviolet renormalon ambiguities in matrix elements don't depend on external states 
and may be calculated at the level of quarks and gluons (Fig. |^). Note that there is 
an ultraviolet renormalon ambiguity in the wave function renormalization AZq 
(Fig. I^d). The result is 



3 Am 
2 m 



km 



2C'a 2 A 
-g^/^mAm, 



Ap: 



3 

mm 



Ap 



1 C'a 2 A 



(45) 



The sum of ultraviolet ambiguities of the 1/m^ contributions to ( pT]) cancels the infrared 
ambiguity of the leading term. 

The requirement of cancellation of renormalon ambiguities in the mass splitting (p8|) for 
all m allows us to establish the structure of the leading infrared renormalon singularity in 
S{u) at M = ^ beyond the large /3i limit. The ultraviolet ambiguity of the square bracket 



m 



should be equal to fi^ times 



A 



V 



me 



27r 



1^2 



a. 



[l + 0(a,)]. 



(46) 



In order to reproduce the correct fractional powers of a^, S{u) in (H0[) should have the 



branch point at m = 1 instead of a pole: 



S{u) 



l+/32/2/3f 



1 19 
2CfKi - -CaK2 + 



CaK. 



12 



-7i/2/3i 



1 

+ 2 



C.K. 



7i/2/3l 



(47) 
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dl 



d2 



Figure 5: Diagrams for ; quark loops are inserted in all possible ways. 



where omitted terms are suppressed as | — u compared to the displayed ones. Normaliza- 
tion constants are known in the large jSi limit only: Ki = 1 -\- 0(1/ j3i). The large-order 
behaviour of the perturbation series for 6c is 

c„+i = n\ {2(3,r n^''^^" {^CpK^ - ^CaK^ + fCAK.n-^^/'^^ + CaK^tC-'^^-] , (48) 

where omitted terms are suppressed as 1/n compared to the displayed ones. 
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